COMMUTING LINEAR OPERATORS AND DECOMPOSITIONS; 
APPLICATIONS TO EINSTEIN MANIFOLDS 



A. ROD COVER AND JOSEF SILHAN 



Abstract. For linear operators which factor P = PqPi ■ ■ - Pi, with suitable 
assumptions concerning commutativity of the factors, we introduce several no- 
tions of a decomposition. When any of these hold then questions of null space 
and range are subordinated to the same questions for the factors, or certain 
compositions thereof. When the operators Pi are polynomial in other commut- 
ing operators then we show that, in a suitable sense, generically factorisations 
algebraically yield decompositions. In the case of operators on a space over 
an algebraically closed field this boils down to elementary algebraic geometry 
arising from the polynomial formula for P. Applied to operators P polynomial 
in single other operator V this shows that the solution space for P decomposes 
directly into a sum of generalised eigenspaces for V. We give universal formulae 
for the projectors administering the decomposition. In the generic setting the 
inhomogenous problems for P reduce to an equivalent inhomogeneous problem 
for an operator Hnear in V. These results are independent of the operator P, 
and so provide a route to progressing such questions when functional calculus 
is unavailable. Related generalising results are obtained as well as a treatment 
for operators on vector spaces over arbitrary fields. We introduce and discuss 
symmetry algebras for such operators. As a motivating example application 
we treat, on Einstein manifolds, the conformal Laplacian operators of Graham- 
Jenne-Mason-SparHng. 



1. Introduction 

A motivating algebraic question is as follows. For V a vector space, V : V ^ V 
an arbitrary linear operator, and P : V — > V a linear operator which is polyno- 
mial in V, then what do we know about the solution space for P in terms of the 
generalised eigenspaces of V7 The question is obviously most interesting when V 
is infinite dimensional. In fact we want to treat this, and related questions, uni- 
formly without using any information about the operator V or the vector space V. 
Obviously any gains in this direction are particularly important in settings where 
functional calculus is unavailable, but they also provide a potentially important 
first simplifying step even when there is access to functional calculus. 

In the case that the field involved is algebraically closed we obtain a complete 
answer to the question above. 

Theorem 1.1. Let V be a vector over an algebraically closed field ¥ . Suppose that 
V is a linear endomorphism on V , and P = P\V] : V V is a linear operator 
polynomial in T>. Then the solution space Vp, for P, admits a canonical and unique 
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direct sum decomposition 

(1) Vp = ©LoVa, , 

where, for each i in the sum, Vx^ is the solution space for (T> + Aj)^' (pi G Z>oy) 
with — Aj e F a multiplicity Pi solution of the polynomial equation P[x\ = 0. The 
projection Projj : Vp — >■ Va^ is given by the universal formula [34\ )- 

The cross reference ( l34l ) refers an explicit formula given in the next section. If 
M G V satisfies 

(2) (V + X)Pu = 

and is non-zero then we shall term u a generalised eigenvector for V correspond- 
ing to the generalised eigenvalue —A. Using this language a partial paraphrasing 
of Theorem 11.11 is that the solution space for P is a direct sum of generalised 
eigenspaces for T>. The Theorem above is an immediate corollary of Theorem 13. 4[ 
for the case that P is given as a fully factored expression, this states the situation 
for V over an arbitrary field. Related eigenspace/eigenspectral results follow, see 
Corollary 13.81 

We may also consider inhomogeneous problems Pu = /. In the case that F is 
algebraically closed then, by rescaling, this boils down to a problem for an operator 
of the form 

(3) Pu := (V + Xoy^'iV + Ai)P^ ■■■iV + XeY'u. 

Theorem 1.2. Let V be a vector space over a field¥. Suppose that V : V ^ V is a 
non-trivial linear endomorphism and consider P : V ^ V as in ^ with the Aj G F 
mutually distinct and for i = 0, - ■ ■ ,i,piE Z>i. Let us fix f E V. There is a 1-1 
relationship between solutions m G V of Pu = f and solutions {uq, ■ ■ ■ , m^) G ©^+^V 
of the problem 

(4) (V + \oY"uo = /,■■■, + \,ru, = f. 

In Theorem 13.51 we give the explicit transformation between the two problems. 

Given a linear operator V : V ^ V, let us write 7l(V) to denote the image of 
V and J\f{V) the kernel of V. We may summarise part of the key information in 
Theorem 11.21 and Theorem 11.11 (or more accurately Theorem 13.41) by the following. 

Corollary 1.3. For P : V ^ V as in ^ we have 

i=i i=e 

n{P) = fl 7^((P + A,fO, Ar(P) = Ar((I? + A,f 0- 

Taking the special case of D being ^ acting on the smooth functions of M the 
above recovers much of the standard theory of constant coefficient linear ordinary 
differential equations. Evidently these aspects generalise to arbitrary operators 
V. In fact the above results are just examples from a very general context (not a 
priori related to polynomials) in which we develop considerable theory as below. 

Section [2] introduces various notions of a decomposition for linear operators 
P that factor P = PqPi ■ ■ - Pe, and where the factors mutually commute. Each 
decomposition is based on some level of invertibility; not invertibility of P, nor 
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the components individually but rather of the system {Pq,Pi,--- ,Pi). This is 
given initially in terms of identities satisfied by "relative inverses", see (JG]) and 
(fTGll . However a more intuitive picture may be obtained by diagrams as follows. 
One may construct a canonical complex from the operators Pq, Pi, - ■ ■ , Pe (see the 
diagrams ( fTTl) and (fT2l)). This is the Koszul (cochain) complex for the system 
{Pq, ■ ■ ■ ,Pi). In each case the required invertibility means that a certain class of 
subcomplexes of this is exact. The latter is described in Section [2^41 More than 
this see Theorem I2.12[ This shows that, remarkably, making only assumptions 
concerning the exactness of certain subcomplexes of the full Koszul complex we 
recover almost the entire information of a class of decompositions. 

The basic Koszul complex which underlies several of our constructions is also a 
central tool in the treatment of certain spectral systems for commuting operators: 
the so-called Taylor spectrum [23l [10], and the related split spectrum. See [T9] . 
and references therein, for further discussion. Operators polynomial in another 
operator, as above, have also been subject of spectral theory, e.g. [18|. While we 
believe there is considerable scope to develop spectral theory based around our 
discussion, this will be deferred to a later treatment. Our current focus is the 
use of purely algebraic considerations which may be applied rather universally. 
In particular at no point do we need a Banach structure on the vector spaces 
or (spaces of) operators involved. We include some minor comments concerning 
eigenspectrum and eigenspaces. 

For any P admitting a decomposition, of some type, and for any inhomoge- 
neous problem, we are able to describe completely and explicitly the structure of 
the solution space in terms of data for the component operators or certain prod- 
ucts thereof. In particular the results above generalise immediately, see Theorem 
12.21 Corollary 12.31 Theorem 12.101 and Corollary 12. Ill It is meaningful to say that 
the required invertibility for the system (Pq? -Pi? ■ ■ ■ ? Pi), in order to obtain some 
decomposition, is very weak (cf. Chapter IV, Theorem 4 [l9]). In fact for linear op- 
erators polynomial in commuting endomorphisms T>q,'Di, - ■ ■ ,Dk, via elementary 
algebraic geometry we show that it is attained generically. In any specific case, 
over an algebraically closed field, establishing any of the class of decompositions 
boils down to verifying that collections of algebraic varieties determined by com- 
binations of the factors have no common point. See Theorem 13.101 For example 
constant coefficient inhomogeneous linear partial differential equations may gener- 
ically be reduced to equivalent lower order equations using Theorem 12.101 and 
specific problems are practically treatable. We should also point out that for oper- 
ators polynomial in commuting operators the decompositions we study and obtain 
are constructed by purely algebraic means. This means the results we obtain are 
universal; they are independent of the operators Vq,Vi,--- ,1)^. The "relative 
inverses" are given by polynomials in the same operators Vo,Vi, ■ ■ ■ ,Vk. (So for 
example if the are differential operators then the entire theory is within the 
category of differential operators polynomial in these.) Regarding the weakness of 
the relative invertibility conditions see also Chapter IV, Theorem 4 of [T9] . 

For operators P polynomial in a single other operator V, as above, a polynomial 
factorisation of P is generically the strongest (non-trivial) form of decomposition. 
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Thus, and for other reasons, the strongest results are obtained in this setting. 
Some are summarised above and Section 13.11 develops the full theory. 

A non-trivial application for some of these ideas is the study of differential op- 
erators polynomial in the Laplacian. Problems of this nature arise in differential 
geometry and, in particular, in the study of conformal Laplacian type operators. 
One of the simplest examples is the conformal Laplacian Y. This curvature modifi- 
cation of the usual Laplacian A is, in a suitable sense, conformally invariant and its 
importance was observed early last century, see e.g. [6]. Paneitz constructed a con- 
formal operator with principal part [20] and then cubic operators are due to T. 
Branson and V. Wiinsch. Later Graham, Jenne, Mason and Sparling (GJMS) p!5] 
extended these results to a very large family which in odd dimensions, for example, 
includes conformal Laplacian operators of all even orders. Recently this family has 
been seen to play a deep role in Riemannian, pseudo-Riemannian and conformal 
geometry. For example the operators have a central role in the geometry of the 
asymptotically hyperbolic Einstein-Poincare metric which underlies the AdS/CFT 
correspondence of physics, see e.g. [UlIIl]. In another direction the GJMS opera- 
tors control the equations for the prescription of Branson's Q-curvature, and also 
the prescription of the non-critical Q-curvatures [HI [7]. These problems generalise 
the celebrated Yamabe problem (see [21] and references therein) of scalar curva- 
ture prescription. It was shown in [12] that on Einstein manifolds the Q curvature 
and the non-critical Q-curvature are constant. In both cases this result is related 
to another result in [12], namely that on Einstein manifolds the GJMS operators 
are given by factored polynomials in the Laplacian. In section [5] we will use the 
Theorems above with this polynomial factorisation to discuss in any signature, 
and on any non-Ricci-flat Einstein manifold, the relationship of the solution space 
the GJMS operators to the spectrum and solution space of the conformal Lapla- 
cian operator. Via Theorem 11.21 the differential order 2k inhomogeneous problem 
PkU = f for these operators may be reduced explicitly to an equivalent second 
order (Laplacian) problem of the form 

(A + X)u = f 

where A := diag(Ai, ■ ■ ■ , \k) (with the Aj given explicitly in terms of the scalar 
curvature), u = Transpose(Mi, ■ ■ ■ , Uk) and / here means Transpose(/, ■ ■ ■ , /), see 
Proposition 15.41 (In fact, as commented after Proposition 15.41 by enlarging the 
space on which operators may act, the problems in most cases may be reduced 
in the same spirit to equivalent first order problems.) This may have applica- 
tions in the understanding of Q-prescription on conformally Einstein manifolds as 
such inhomogeneous problems may be viewed as model linear problems for the 
true (non-linear) prescription problems. In fact the prescription problems involve 
equations of the form PkU = h{u)f, for h a suitable function of u (in many cases 
simply h{u) = constant. for a suitable power p). The tools of Theorem 13.51 still 
apply when we have a non-linearity of this type and so such equations reduce to 
(A + A)^* = h{bu)f, on non-Ricci-flat Einstein manifolds, where 6 is a row matrix 
of constants determined by the scalar curvature. In yet another direction the tools 
of section 13.11 show that on Einstein manifolds, of any signature, the eigenspaces 
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and eigenvalues of the Pk all arise from generalised eigenspaces and generalised 
eigenvalues of Y. This is just a special case of the general result in Corollary 13.81 
and should have application in the representation theory of the orthogonal groups 
SO{p + 1, g + 1) (which, modulo issues of covering, act as the conformal group on 
appropriate products of spheres) . The case of conformal Laplacians as discussed 
here is just an example application. There are many other settings where these 
ideas apply. For example parallel to the theory of conformal Laplacians there is a 
theory of sub-Laplacians in CR geometry [I3]. Product manifolds yield commuting 
operators (such as Laplacians of the components) and so the machinery of Section 
13.31 is ready for these. 

Recently there has been a growth in interest in the very old problem concern- 
ing the so-called symmetries and symmetry algebras of Laplacian type operators, 
see [HI |9] and references therein. Such symmetry operators play a central role in 
separation of variables techniques for the solution of the Laplacian operators in- 
volved. In section [4] we introduce symmetry algebras which generalise this notion 
to a large class of linear operators. Using this we obtain, for example, general 
results relating the symmetry algebra for a linear endomorphism V and that of 
a second operator P polynomial in V. See in particular Theorem 14. 1[ In section 
15.11 Einstein manifolds are once again used to illustrate these ideas in a concrete 
setting. 

Finally we point out that the theory of polynomial operators may also be applied 
to large classes of differential operators which are not simply polynomial in another 
operator V. This is achieved by, for example, composing operators which do have 
the form P[D] (i.e. P polynomial in a differential operator T>) with other suitable 
differential operators. For example large classes of conformally invariant operators 
on tensor and spinor fields arise this way [4l[22]. This will be taken up elsewhere. 

The authors are grateful to John Butcher, Andreas Cap, Mike Eastwood, V. 
Mathai, Paul-Audi Nagy and Jan Slovak for helpful discussions. The first author 
would like to thank the Royal Society of New Zealand for support via Marsden 
Grant no. 06-UOA-029. The second author was supported from the Basic Research 
Center no. LC505 (Eduard Cech Center for Algebra and Geometry) of Ministry of 
Education, Youth and Sport of Czech Republic. 

2. The general setup 

Our study here will concern linear operators given by a composition P = 
PqPi ■ ■ ■ Pe where the factors mutually commute. In the case that each factor is in- 
vertible then the essential properties of P are given simply in terms of the factors. 
Otherwise the situation is significantly more complicated. Here we explore part of 
the latter domain. In general there are relative qualities of linear operators Pq cind 
Pi that dramatically affect the nature of the composition PqPi- As a very simple 
example, and a case in point, one may compare (^ + A)(^ + /i), with A 7^ /i, 
to (^ + A)^. As operators on the line these are rather different beasts. These 
observations in part motivate considering the following class of linear operators. 
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2.1. Decompositions of linear operators. Let V denote a vector space over a 
field F. Suppose that P : V ^ V is a linear operator which may be expressed as a 
composition 

(5) P = PoP,---Pi 

where the linear operators Pi : V ^ V, i = 0, ■ ■ ■ , i, have the following properties: 
there exist linear operators Qi: V—>-V,i = 0,l,---,i, that give a decomposition 
of the identity, 

(6) idv = QoP° + ■ ■ ■ + QiP' , 

where P' := n^^^^^gPj, i = 0, ■ ■ ■ , i; and the PjS and the QjS are mutually com- 
muting to the extent 

(7) P,P,=P,Pi, and P,Q,=QjP,, ^, j G {0, ■■■,£} ; 

When i = this may be viewed to hold trivially. For other cases we shall use the 
following terminology. 

Definition. For a linear operator P : V — > V an expression of the form ([5]) will be 
said to be a decomposition of P if the factors Pi, i = 0, ■ ■ ■ ,i, satisfy the conditions 
just described and i > 1. 

Note that if one of the factors P* is invertible (from both sides) then we have 
([6]) immediately. In general requiring the identity ([6]) is a significantly weaker 
requirement. It states for example that the operator (P°, ■ ■ ■ , P^) : ®1=qV — > V has 
a right inverse given by the operator {Qq, • • • , Qi) : V — > Q)1=qV. (This statement 
also holds if we swap the roles of the PjS and the QiS.) 

We next observe that the identity ([6]) controls a decomposition of the null space 
AA(P). 

Lemma 2.1. For each i E {0, • ■ ■ ,£}, we have 

Q,P':N{P)^N{P,) 

and this is a projection. 

Proof. Since PiQiP' = QiP it is clear that QiP' : ^^{P) ^^iP^). Then on V, 
and hence in particular on Af{P), we have the identity given by (l6|). But for j ^ i, 
Pi is a factor of P^ and hence P^ annihilates N{Pi). So QjP* is the identity on 

M{Pi). □ 

For convenience we will often use Vp to denote the null space of a linear operator 
P on V, so e.g. we may write QjP* : Vp — > Vp. . 

We consider now the inhomogeneous problem Pu = f. Of course the solution 
space is the affine subspace in V obtained by translating Vp (the solution space 
for the linear problem) by any single "particular" solution to Pu = /. It turns out 
that, by applying ((61) to HQTl^Pi), we can decompose the inhomogeneous problem 
to a simpler inhomogeneous problem in a way that generalises the treatment of 
the homogeneous cases. 
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Theorem 2.2. Let V be a vector space over a field ¥ and consider P : V V as 
in ^ with the factorisation there giving a decomposition, i.e. ^ and ^ hold. 
Let us fix f G V. There is a 1-1 relationship between solutions m G V of Pu = f 
and solutions {uq, ■ ■ ■ , Ue) G ©^^^V of the problem 

(8) p^uo = fr-- ,PeUe = f. 

Writing Vp for the solution space of Pu = f and (for i = 0, ■ ■ ■ , i) for the 
solution space of PiU = f. The map F : Vp ^ x^^qV/ is given by 

u ^ {P\, ■ ■ ■ , P^u) , 

with inverse B : x^^qV/ — > Vp given by 

1=1 

{uo, ■ ■ ■ ,ue) h^^Q 



Ui . 



i=0 



On V we have B o F = idy, while on the affine space x^^qV/ we have F o B = 

id e yf. 

Proof. Suppose Pu = f. Then PiP^u = Pu = f and so Fu is a solution of (I8|). For 
the converse suppose that {uq, ■ ■ ■ , m^) is a solution of ([8]) and write u := X]i=o 
Then 

i=£ i=e i=e 

Pu = Y, PQi^^ = Yl Q^P'P^^^ = H Q^P'f = f 

i=0 1=0 i=0 

where finally we have used ([6]). 

It remains to establish the final assertion. By construction B o F = ^^Iq QiP^ 
and so BoF = id^ is just the identity ([6]). (Then in particular BoF = idys .) Next 

we calculate Fo 5 on xf^gV/. For the /c*'^-component, we have [FB{uq., ■ ■ ■ ,ui)\k 
given by 

i=i 



i=0 

Using the commutativity of terms, and that PiUi = /, this gives 

i=£ 171=1 
( E n Prnf)+QkP' 



kj^i=0 i,kj^m=0 

Now using ([6]) and then PkUk = /, we obtain for the last term, 

i=£ m=e 

QkP^Uk = Uk- Y Qi W Pmf ■ 

ky^i=0 i,kj^m=0 

Thus 

[FB{uo, ■ ■ ■ ,ue)]k = Uk , 
for any G {0, ■ ■ ■ ,i} and we conclude that FB is the identity on x^^qV/. □ 
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For any operator of the form ([5]), with the Pi mutually commuting, we obviously 
have +^=^7V(Pi) C A/'(P) and 7^(P) C ^^=^7^(Pi). From the above we see that 
if (l6|) holds then these containments are equalities. In summary we have the 
following. 



Corollary 2.3. For P : V ^ V, with ^ giving a decomposition, we h 



ave 



np) = f]np,), Ar(p) = 0Ar(Po. 

i=0 i=0 

The decomposition of N{P) is given by the identity 



i=0 

where, for each i E {0, - ■ ■ Projj : Af{P) J\f{Pi) is the projection given by 
the restriction of QiP^ from 

In section im we shall show that operators polynomial in a single other operator 
generically admit decompositions that may obtained algebraically and explicitly. 
The explicit formulae for the Q's (in the identity (l6|)) are given in terms of the 
basic data of the factorisation PqPi ■ ■ ■ Pi- Applications treated in Sections [4] and 
[5] then show that the decompositions are a powerful tool. Before we take these 
directions we study the algebraic structures underlying a decomposition and this 
leads to results which show that the decompositions are a special case of a rather 
general theory with a vastly broader scope for development and applications. 

2.2. Relative invertibility, and operator resolutions. We first shed some 
light on the above constructions. 

In relation to the identity ([6]), suppose that we have linear endomorphisms 
Po, Pi, on a vector space V, and there exist further endomorphisms Qi, Qo so that 

idv = QiPo + QoPi. 

Then clearly Pq is injective on the null space of Pi and is an invertible endomor- 
phism of A/'(Pi) space if, for example, the commutativity relations ([7]) hold. So 
the displayed identity manifests what we might call relative invertibility of the 
operators Pq and Pi. Evidently we may solve Pqu = / for / G Af{Pi). This is a 
consequence of the fact that the short complex associated to the system, 

^'0 V 

(9) V ^ © V 

V 

is forced to be exact (and is split) by the identity id\; = QiPq + QqPi. The splitting 
sequence takes the same form with Qo and —Qi formally replacing, respectively, 
Pq and Pi. For example the system 

PqU = /o, PiU = fi 
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has the exact integrability condition Pi/o = Pq/i and if this holds then the solution 
u = Qifo + Qofi is unique. Including the projections for the bundle in the centre 
of the sequence ([9]) we obtain a diagram 



I V I 

— ^^C^ — 

Here the long arrow indicates the composition PiPq and note that by viewing the 
centre column as a direct sum we include the information of the original complex. 

A related observation follows. This concerns how, for the very simple case of 
£ = 1, the essential content of Theorem 12.21 is captured in the short exact sequence 
([9]). Notation is as above. 

Lemma 2.4. // ^ is exact and (m^, u^) solves the system PqvP = f and PiU^ = f, 
then {u^, u^) = {PqU, Piu) for a unique m G V satisfying Pu = f . 



Proof. Note that is in the null space of (Pi — Pq) and so, using that ((91) 

is exact, we have the result. □ 



A key point is that this holds without explicit mention of the splitting Q-operators. 
Of course, for example, (Pi — Pq) has a left inverse but we do not any commuta- 
tivity properties of this beyond what is forced by ([9]) being exact. We will return 
to this point in Section 12.41 
The case i = 2, i.e. the system 



(10) 



PqU = /o, Piu = /i, P2U = /s 



demonstrates the general situation more accurately. Here we have P = P0P1P2 
where the commutators [Pi,Pj] are all trivial. If there is a solution to (fTOjl then 
it is necessary that Pifj = Pjfi for all < i, j < 2. These and further problems 
with their integrability conditions may be organised into the complex 

(11) 




Now consider for each operator — ^ also a corresponding operator < — - in the op- 
posite direction. We assume, as before, [Qi,Pj] = for i,j G {0,1,2}. (Note 
that if also the operators Qi are mutually commuting then they too form a com- 
plex.) Denoting the space of the degree p G {0,1,2,3} by V(p) (i.e. V(0) = V, 
V(l) = V © V © V etc.), these complexes become 







V(0) 



Q(i) 



V(l) 



0(2) 



V(2) 



Q(3) 



V(3) 
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where the operators P{p) : V{p) — * V(p + 1) are given by the corresponding sum of 
operators ±Pj (and similarly for Q{p)). Now we can write the system (fTOll simply 
as P{0)u = f where / = (/o, f\, /a) G V(l). 

The cohomology of this complex is related to the solution spaces of the problems 
P{p)u = 0. In particular = J\f{P{0)). While in general little could be said 
about the cohomology the key point is this. Observe that 

Q{p + l)P{p) + P{p - l)Q{p) = QoPo + Q,P, + Q2P2 

on V{p) for p e {0,1,2,3} (where P(-l), P(3), Q(0) and Q(4) are indicated 
trivial mappings). Hence if the right hand side of the last display is equal to 
the identity then P{p — l)Q{p) is the identity on M{P{p)) and so the complex is 
exact; in the case that the P(p)-complex is exact we shall say the complex is a 
resolution of the operator P. When we have such a resolution then, for example, 
the problem P(0)m = /, has a solution only if we have the "integrability condition" 
/ G N{P{1)) and if this holds then the solution is unique. 

Note that the diagram (fTTl) is essentially the Hasse diagram (or lattice diagram) 
for the natural poset structure of the power set 2^ in the case L = {0, 1, 2}. The 
situation for a general £ G N is analogous, and we shall exploit the connection to 
poset structure to organise the notation. So we consider operators Pj : V — *• V, 
i G L := {0, !,...,£} which are mutually commuting, but otherwise arbitrary. 
The complex will be constructed using 21^' copies of V as follows. The copies of V 
will be indexed by subsets J C L, i.e. Vj := V, and we define the operators 

Pj,i ■■= (-1)I-'<'IP. : Vj ^ Vjuw, JCL,teL\J 

where 

\J <i\:=\{j e J \j <t}\, JCL. 

Further we put 

Vip) := Vj and P(p) := +,=|^|,^^P,, : V(p) ^ V(p + 1). 

JCL, 
\J\=p 

Proposition 2.5. The operators P{p), < p < i form a complex 
(12) — V(0)-^ •■■ -^V(f+1) — 0. 

Moreover, if idy = QqPq + . . . + QePe for some operators Qi : V V , satisfying 
[Qi, Pj] = for i,j = 0,...,£, then this complex is exact. 

If the complex (fT2ll is exact we shall call it a resolution of the operator P. In the 
treatment of the Taylor spectrum for commuting operators on a Banach space the 
Koszul complex here is said to be Taylor regular if it is exact. The main part of the 
Proposition here is in Proposition 3, Chapter IV of [l9|. We include the proof here 
to keep the treatment self-contained and in terms of a single notational system. 

Proof. First we need to show that P{p + 1) o P(p) = 0. This map is a sum of 
mappings Pj,i,j : Vj — > Vj\j{ij} such that \J\ = p and i,j G L \ J given by the 
restriction of P{p + 1) o P{p) to the source subspace Vj C V{p) and projection 
onto the target subspace Vju{i,j} ^ V(p + 2). Fix such a triple {J,i,j) and assume 
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i < j. Denoting qi = \ J < i\ and = | J < j|, we obtain that Pj^ij is the sum of 
the two composite operators Vj — > Vju{i,j} in the following diamond: 

(13) Vju{i} 




But from this we see immediately that Pj^ij = 0. 

Now assume we have operators Qj : V ^ V so that idy = QqPq + • • • + QiPe 
and [Qi, Pj] = as in the Proposition. Consider for every operator — ^ also the 

operator in the opposite direction. (Then the operators labelled by ±Qi also 
form a complex, provided [Qi, Qj] — but we will not need this fact.) We obtain 
the diagram 

P(0) P(l) ^'(^-i) 

V(0) V(l) ^^V(£)^^V(£+1)^^0 

Q(l) Q{2) Q{e) Q{e+1) 

where Q{j) is the sum of the between the corresponding subspaces of V(j) and 
V(j-l) for j = + We denote by P(-l), P{£ + 1), Q{0) andQ(£ + 2), in 

an obvious way, the trivial operators at the left and right extremes of the diagram. 
Let us fix p G {0, ...,£ + 1} and consider the restriction Q{p + 1) o P(p)|vj for 
some J <Z L with \J\ = p. By definition, P{p)\Vj is the sum of operators Pj^i for 
i ^ J. The Q-operators from Vju{i} ^ V(p+ 1) (i.e. the target space of Pj^i) back 
to V{p) correspond to j e JU{i} and have V(ju{j})\{j} Q V(p) as the target space. 
For a given i ^ J, the choice j := i yields the composition QiPi : Vj — > Vj, and 
the choices j & J yield the operators 

(_l)|J<i|p. (_i)IUu{i})\0}<i|Q 

(14) Rij = Vj ^ Vju{i} ^ ^{Ju{i})\{j} 

= qijPiQj ■ Vj V(ju{i})\{j}, J,j e J. 

where Qij e {+1; ~1} is determined by the previous display. Summarising, we 
have obtained 

g(p + i)oP(p)|v, = (5^g,p,) + Yl R^3■ 

\ / on Vj 

igj i^J,j£j 

The same analysis of P{p — 1) oQip) yields 

p(p-i)og(p)|v,= (J^p.g.) + 5^ i?: 

^ . _ ^ on v_i , _ 



where 



{-l)U\{3)<3\Q. (_l)IAO}<i|p. 

^ ^J\{3} 



'15) ^ij^ ^^A{j} ^^{A{j}Mi} 
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Summarising again (and using [Pj, Qj] = 0), we obtain 

Q{p + 1) o P{p) + P{p - 1) o Q{p) Iv, = 



onvj 

i=0 i<^JJ<^J 
t 

The first sum is the identity according to the assumption. To compute the second 
one we use the explicit form of qij and q'^j given by respectively (fT4l ) and f fTSD . If 
2 > j then 

(^_iyA{i}<il = (-_i)l(^u{i})\{i}<j| (•_iy^<*l = _(^_i)IA{i}<'l 

hence qij = —q'ij- One easily sees the latter is true also for i < j. Therefore we 
obtain 

Qip + 1) o P{p) + P{p - 1) o Q{p) = idv(p), 
whence P(p — 1) oQ(p) is the identity on Af{P{p)) and the Proposition follows. □ 

Note that the identity idy = Yll=o QiPi is in general far weaker than ([6]) required 
for a decomposition. This motivates a rather broader notion of decomposition that 
we now introduce. 



2.3. General case: a-decompositions. We define and discuss here a generali- 
sation of the notion of a decomposition which has the decomposition from Section 
12.11 as simply an extreme (but important) special class. Consider the operator 
P = Pq - ■ ■ Pe from (I5|) and the power set 2^ of the index set L := {0,1, ... 
We shall use the notation Pj := YljeJ for 7^ J C L and set P$ := idy. Now 
choose a nonempty subset a C 2^ and assume there exist operators Qj : V — > V, 
J E a that give a decomposition of the identity 

(16) tdy = Y,QjP-^ 

Jece 

where P'^ = Pl\j, and PjS and Qjs satisfy 

(17) PiPj = PjP, and PiQj = QjPi ieL,Jea. 

Definition. For a linear operator P : V ^ V, an expression of the form ((51) will 
be said to be a a -decomposition of P if the identity ( flGl l holds with (fTTll satisfied 
and ^ a C 2^, L ^ a. 

The case of (I5l) being a decomposition is a special case of an a-decomposition 
with a = {JCL||J| = l}. Toward understanding a-decompositions we employ 
a dual notion of a decomposition, as follows. 
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Definition. We say that P = Pq - ■■ Pi is the dual P -decomposition, 7^ /? C 2^, 
{0} 7^ /? if for every J G /? there exist operators Qjj G End(V), j E J such that 

(18) idv = 5^ Qj,,P„ [P„Pfc] = [Qj,„P.]=0, ^,A;gL,jgJ. 

Each system a C 2^ is partially ordered be restricting the poset structure of 
2^. The sets of minimal and maximal elements in a will be denoted by Min(a) 
and Max(a), respectively. We say the system /5 C 2^ is a lower set, if it is closed 
under taking a subset. (That is, if / G /3 and J ^ I then J G 13.) The upper set 
is defined dually. The lower set and upper set generated by a system a C 2^ will 
be denoted by C{a) := {J C / | / g a} and lA{a) := { J ^ / | J ^ L and / G a}, 
respectively. 

Lemma 2.6. Let a C 2^. Then P = Pq - ■ ■ Pg satisfies the following: 

(i) it is an a -decomposition <^==^ it is a Max{a) -decomposition <^==^ it is an £(«)- 

decomposition 

(a) it is a dual a -decomposition it is a dual Mm{a) -decomposition <^=^ it is 
a dual U (a) -decomposition. 

Proof. The proof of (i) follows easily from the definitions, the proof of (ii) is also 
obvious. □ 

To formulate the relation between a- and dual a-decompositions, we need the 
following notation. We put a" := 2^ \ C{a) and := 2^ \ U{a). Clearly [d^Y = 
C{a) and (a')" = lA{a). Also it is easily seen that 

a" = {J C L I V/ G « : J\ / ^ 0} 

^ ' = {J CL\WI ea: I\J y^ij)}. 

The first part of the following proposition describes the duality in the special 
case dll). 

Proposition 2.7. ^ is equivalent to 

(20) idv = Qi,jPi + Qj,iPj 

where Qi,j G End(V) and satisfy [Qij,Pk] = for every triple of integers {i,j,k) 
such that < i, j, k < i and i j . That is, with is equivalent to the dual 
P -decomposition for P = {J C L \ \ J\ = 2}. 

More generally, ^ is an a -decomposition if and only if it is a dual a"^ -de- 
composition. Equivalently, ^ is a dual P -decomposition if and only if it is a 
P^ -decomposition. 

Proof. We shall prove the first part of the general statement, i.e. that (l5|) is an 
a-decomposition if and only if it is a dual ^"-decomposition. Also we will suppose 
a = C{a). This is no loss of generality due to Lemma [2?6l 

Assume ([5]) is a-decomposition and consider J G a". That is, J \ / 7^ for all 
/ G a. For any J C L, we have L \ I = J' U {J \ I), J' C L as the disjoint union. 
Hence 

QjP' = QiPl\i = {QiPj')Pj\i. 
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From this, it is obvious that the identity idy = ^/^o Qi^^ be easily rewritten 
to the form (fTSll because 7^ J \ / C J for all J G a. (The commutation relations 
in (fTSll are clearly satisfied.) 

Now assume ^ is the dual a"-decomposition, i.e. idy = Ylij<^jQjjPj fo'^ every 
J G a". We shall prove that for every J G a" we have a decomposition of the 
identity 

(21) idv = J2 Qi,jPj\i^ Qi,J e End(V), 

such that [Qi,j, Pk] = for every / G a, and k E L. Then the proposition follows 
from the choice J := L G a" as Pl\i = ■ The proof will use induction on the 
partial ordering of 2^ (given by inclusion), will use that 2^ is the disjoint union 
2^ = a" U a, and also that, since a = C{a), we have the least element G a. 

Before we do the induction let us first consider as easy case which indicates how 
the argument works, viz. J G Min(a;"). It follows from this minimality that for 
every j G J we obtain Ij := J \ {j} G a hence {j} = J \ Ij for some Ij G a. 
Using this and since (l5|) is a dual a"-decomposition and J G a", we conclude 
idv = Y.jejQjjPj = '^j€jQjJ^J\ij- t^*^ latter sum is of the form (ET]) 

because Ij G a. (We put Qi^j := for every J G a not of the form Ij for some 
j G J.) The commutativity conditions in (i2Ti ) follow from the definition of the 
dual a"-decomposition. 

Now consider J G a". Since 2^ = a^Ua is a disjoint union, there are sets J' and 
J" so that J = J'U J" where J \ {j} =: Ij G a for j G J', and J \ {j} =: Jj G a"" 
for j G J" . Now, as J G a", the assumption of a dual a"-decomposition gives the 
identity 

idv = J2 Qj,jPj\i, + Yl '^jj^i , 

jeJ' j<^J" 

where we have used that J \ Ij = {j} for j G J'. This first sum is of the form 
required in (l2Tl l. the second one is not. But since a" 3 Jj <Z J for j G J", we may 
assume, by the induction, that idy = ^leaQjjJ^JjV- Acting on this by Pj, we 
obtain Pj = '^j^^ Qjj^i^AiiMj})- Here / \ {j} G a (because a = C{a)) hence the 
latter sum is of the form on the right hand side of (i2Ti ). Consequently, putting 
these expressions for Pj, j G J" into the previous display, we obtain decomposition 
of identity of the form of (|2T]1 . The required commutativity relations clearly hold 
thus the proposition follows. □ 

Remark 2.8. A main point of the Proposition above is to shed light on the nature 
of a-decompositions. The a-decomposition is what gets directly used in studying 
the solution space for P. However at first this seems rather mysterious since, for 
example, the P'^ in the identity (fTBl l are complementary to the Pj. The first part 
of the Proposition exposes one view of what it means to say that PqPi • • ■ is a 
decomposition: it shows that ([6]) is equivalent to the PjS being mutually relatively 
invertible. We will see in Section [2^41 that this picture generalises. 

Next note that the proof above, begining with (ITSll . inductively constructs ex- 
plicit formulae for the Qj in (fTGll in terms of products of the Qjj from (fTSll . Note 
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also that although, we do not require [Qjj, Qj',j'] = for J, J' G a, j G J, j' G J' 
in (fTSll . in the special case (l20ll in the Proposition one shows that [Qi,j,Qj,i] = 
easily follows from (l20l) and the vanishing of the [(5ij,-Pfc] as assumed. 

The subsets a C 2^ are partially ordered by inclusion (i.e. now we use the poset 
structure of 2^ ). Given an operator P in the form ([5]) consider the family F 
of systems a such that (I5l) is a dual a-decomposition. Then T has the greatest 
element ap = IJaer*^- Then an "optimal" choice for the (dual) a-decomposition 
of P is a := Min(ap). (We want to have in a to the smallest possible subsets of 
L. So if the PjS are not invertible then the case of a dual decomposition may be 
regarded as the best we can do. With this philosophy we thus take ap. Then using 
Lemma [2.61 we take a := Min(ap) as it is easier to work with a smaller number 
of subsets.) Consequently, we obtain the optimal choice /3 := Max((ap)') for the 
/^-decomposition of P. 

In the case one is able to decide, given a subset J L, whether idy = 
Ylij£jQj,i^3 for some Qjj G End(V), it is easy to find the optimal (dual) de- 
compositions. This is, for example, the case of polynomial operators discussed in 
Section 13.31 

The decomposition used in Lemma ETT] and Theorem 12.21 is a special case of the 
a-decomposition, 7^ a C 2^ where L = {0, 1, . . . ,£}. The null spaces of the Pj 
will in general meet non-trivially. However note the following. 

Lemma 2.9. //a C 2^ gives an a -decomposition of P then 

(22) QiP^ : 7V(P) ^ 7V(P/) for all I e a . 

If a satisfies I H J = (I) for all I J E a then, for each I E a, QiP^ in [2^) is a 
projection. 

Proof. The point is that if the sets in a are mutually disjoint then P"^ (and hence 
QjP'^) annihilates Af{Pi) whenever I ^ J . So the proof of Lemma [2TT] generalises 
easily. □ 

Using the Lemma and by an easy adaption of the proof of Theorem 12.21 we 
obtain the following. 

Theorem 2.10. Assume P : V V as in ^ is an a -decomposition. Let us fix 
/ G V. There is a surjective mapping B from the space of solutions {uj)j<za G 
©'"'V of the problem 

(23) Pjuj = /, Jea. 

onto the space of solutions u G V of Pu = f . 

Writing Vp for the solution space of Pu = f and (for J & a) Vj for the solution 
space of Pju = f. The map B : Xj^^Vj Vp is given by 

A right inverse for this is F : Vp ^ x JeaVj given (component-wise) by 

u I— > P'^u ; 
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on V we have B o F = idy. 

If a satisfies I H J = % for all I ^ J ^ a then, F is a 1-1 mapping and F o B 
is the identity on the solution space to ( f^j . 

Hence the generalisation of Corollary 12.31 is as follows. 

Corollary 2.11. For P : V —>■ V, with ^ giving an a -decomposition, 7^ a C 2^, 
we have 

n{p) = fi n{p.j), U{P) = +jeaM{Pj). 

If a consists of mutually disjoint sets then we have 
and this is given by 

where, for each I E a, Projj : Af{P) J\f{Pi) is the projection given by the 
restriction of QjP^ from [TW . 

So although the assumption of an a-decomposition, for an operator P, is in general 
a vastly weaker requirement than that of a decomposition, we still have the critical 
result that one may solve the inhomogeneous problem Pu = / by treating a "lower 
order" problem involving the same inhomogeneous term /. 

2.4. a-decompositions in terms of operator resolutions. Recall that the 
complex (fT2l) in Proposition 12.51 promotes to being an operator resolution (i.e. 
is exact) if we make the the rather weak assumption idy = Ylo Qi^i (with the 
usual commutativity of operators assumed). On the other hand expression (l20l) 
in Proposition 12.71 shows when one has a decomposition P = PqPi ■ ■ ■ Pi (i.e. ((61) 
holds) then every diamond subcomplex (fT3ll of the operator resolution diagram 
(fT2ll is exact. In a sense, that we now make precise, this is the key algebraic 
content of a decomposition. 

Consider then P = PqPi ■ ■ ■ Pi, where as usual the Pi G EndV are mutually 
commuting. We have the complex f |T2i ). Let us assume that in this each diamond 
subcomplex of the form (fT3ll is exact (in the sense of ([9])). Then the complex f |T2i ) 
is exact and so gives an operator resolution. We shall investigate to what extent 
the results for decompositions survive if we take this setting without explicitly 
requiring the identity (l6|). 

We earlier discussed the case i = 1. To shed light on the general situation we 
look now at the case i = 2, P = P0-P1-P2, so we have the complex (fTTl) . For / G V, 
consider the inhomogeneous problem PqVq = /, PiVi = f , P2V2 = f ■ Since each 
diamond is exact we have that 

( Vq \ ^( -P1U2 \ f f P2Uo\ 

\vi J \ -P0U2 J \V2 J \PlUo J - 
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These are consistent only if P2U0 + P0U2 = 0, and, when this holds, using that the 
P07 P2 diamond is exact we find that, for i = 0, 1, 2, f j = for m G V satisfying 
Pu = f (cf. Theorem EIID. 

The results for £ = 1,2 extend to general £ G M. 

Theorem 2.12. Suppose that we have P = PqPi ■ ■ ■ Pe, as in Suppose also 
that in the corresponding sequence flE) every diamond [T3\) is exact in the sense 
of Then all results of Theorem \2.2\ hold except the map B should replaced by 
the map B' given [u^, ■ ■ ■ , u^) u hy taking in ( fl^j the unique preimage (of the 
map F) inVij) of 

- e Vl\{o} © ■ ■ ■ © Vl\{£} 

solving 

(24) Piu' = f, 2 = 0,1, 

We have 

i=i i=e 

n{P) = []n{p,), Ar(p)-0Ar(p,). 

i=0 i=0 

Proof. First note that since each diamond in the sequence (fT2ll is a complex then 
PiPj = PjPi for i,j G {0, ■ ■ ■£} and the sequence is a complex. 

If / G 7^(P) and Pu = f then recall that ■■■ ,u^) := (P°m, ■ ■ ■ ,P^u) is 
a solution of (l2ll) . We will (strong induction to) prove that any solution of (l2ll) 
has this form, as forced by the consistency of exact diagram (fT2l) . Note that by 
Lemma Em this is true for case £ = 1. 

Assume now that i >2. Starting at V{o} and V{i), in the (length resolution 
diagram f |T2ll for Pq ■ ■ ■ Pe, there are subcomplexes of length i that each take the 
form of (fT2ll : in both of these the terminal space is Vl (where, as usual, L := 
{0, 1, ■ ■ ■ By the inductive hypothesis, consistency of these subcomplexes 

mean that there is Uq G V{o} satisfying P°Mo = / and similarly P^Ui = f. (Recall 
P* means P/Pi ) Now V{o,i} is in both subcomplexes and we obtain a consistency 
condition: by the process of repeatedly using Lemma 12.41 to take preimages and 
enforce consistency (at each diamond) in order to solve for uq and ui, it follows 
easily that Uq and Ui are both "potentials" for the (by induction unique) entry in 
V{o,i}. Hence 

(-^1 Po)(;;°)=OGV,o,). 

Since the diamond (fT3l) for J = is exact (in the sense of ([9])) it follows that 
necessarily Uq = PqU and Ui = PiU for some u G V0. From these it follows, 
respectively, that = P*m, for i = 1, ■ ■ ■ , £ and = P^u, for j = 0,2,3, ■■■ ,£. It 
also follows that Pu = f. 

By construction B' is 1-1 and F o B' is the identity on the solution space. That 
the forward map (in the notation of Theorem 12. 2p P is 1-1 is an easy consequence 
of the injectivity of {Pi,Pj) : V ©^V for each pair distinct pair {i,j) & L x L. 
The direct sum in last display follows for the same reason. □ 
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It seems likely that there are analogous simplifications for the general a-de- 
compositions. It has also not escaped our attention that these ideas suggest that 
there should be extensions of the ideas here to the setting where one has a suit- 
able commuting diagram but without assuming that diagram is constructed from 
commuting operators. This will be taken up elsewhere. Note that although the 
Theorem here is conceptually powerful and a far stronger result overall than The- 
orem [221 it seems likely that in practice the identity ([6]) with (I7|) is rather useful. 
In particular one then obtains the projections in Lemma 12. 1[ Also, as we shall 
see in the following sections, for a large class of operators we have have all these 
identities algebraically. 

The resolution diagrams give us a "pictorial" understanding of the a-decomposit- 
ions. For each dual /3-decomposition P = PqPi ■ ■ ■ Pe and J E (3 we have idy = 
Yljej QjjPj (with appropriate commutativity conditions) and so a collection of 
length I J| exact subcomplexes of the resolution for P. Each of these is itself an 
operator resolution for HjeJ-^i- ^'^^ ^^^^ of the |J| as we range over J G /5 gives 
some measure of the strength of the dual /^-decomposition: the smaller the sets 
J E a the stronger the decomposition. For example P = {{0}, {1}, ■ ■ ■ {£}} is 
the case that all the Pi are invertible. The duality in Proposition 12.71 allows us 
therefore to understand a-decompositions in the same way: Small sets J in a 
indicate a strong decomposition. 

Remark 2.13. Note that the complexes f |T2D discussed in Section [2^2] were con- 
structed from an arbitrary set Pq, . . . ,Pe of mutually commuting endomorphisms 
of V. Hence using the notation used in (flGl) . we can take this set to be {P'^\J G a} 
for some nonempty system a G 2^. Then it follows immediately from the above 
proposition that if ((51) is an a-decomposition then the corresponding complex is 
exact. 

Remark 2.14. As a final point we note that there are other approaches to the 
inhomogeneous case that naively seem similar to Theorem 12. 2[ For example note 
the following. Assume P to be in the form ([5]) (with the factors not necessarily 
commuting) . Then clearly Pu = f has a solution if and only if there is a sequence 
/o, • • • , e V satisfying 

(25) Po/o = /, Pifi = /o, . . . , Fife = h-i- 

So it is sufficient to find such a sequence to obtain a solution u = fe of Pu = f. 
However this is simply a variant of the idea from differential equation theory where, 
through the introduction of new variables, one replaces a differential equation by 
a system of lower order equations. This is very different from Theorem 12. 2[ The 
system here does not replace Pu = f with a new inhomogeneous equation, but 
rather replaces it with a sequence of problems. We do not have the "source term" 
/o in Pi/i = /o until we have solved the previous problem Pq/o = / and so on. 

3. Algebraic decompositions 

Here we consider operators P polynomial in mutually commuting operators 
Dq, ■ ■ ■ , Dfc. In this setting we show that generically we obtain a-decompositions. 
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In fact in this Section we derive those decompositions (and a-decompositions) 
that may be obtained in a purely algebraic or algebraic-geometric manner from 
the polynomial formula for the operator. Thus these are universal results that are 
independent of the operators Dq,--- , D^. An important feature of these cases is 
that the "relative inverses", viz. the Q-operators in (l6j) and Theorem 12. 101 are then 
also obtained as operators polynomial in the same operators Dq, - ■ ■ ,Dk. Thus if, 
for example, we dealing with P a differential operator then these relative inverses, 
are also differential operators. 

The simplest setting and the strongest results are obtained in the case of oper- 
ators polynomial in a single other operator. Here we derive explicit formulae for 
the decomposition that are significantly simpler and more efficient than expected 
from the general setup. 

3.1. Operators polynomial in a single operator V. Let V be a vector space 
over the field F. Suppose that V : V ^ V is a non-trivial linear endomorphism. We 
may consider the commutative algebra ¥[D] of consisting of those endomorphisms 
V — > V which may be given by expressions polynomial (with coefficients in F) in 
T>. Clearly there is an algebra epimorphism from F[a;] onto ¥[D] given by mapping 
a polynomial P[x] = X^iLo '^i-'^^ operator Pi'D], a formula for which is given 

by formally replacing the indeterminate x in P[x] by T>. That is, a formula for 
P[D] : V ^ V is X]i=o where we write as a shorthand for the i-fold 
composition of V. This algebra map sends 1 G F[a;] to id^. 

We begin by treating operators of the form ([3]). That is P = PqPi ■ ■ ■ Pe where 
Pi = {V + Ai)^% with the Aj G F mutually distinct and for z = 0, ■ ■ ■ , ^, pi G Z>i. 
Since the algebra ¥[D] is commutative, we may access the results of Section 12.11 
provided we obtain the identity ([6]). This we have from the Euclidean algorithm 
as follows. To a polynomial of the form 

(26) P[x] = ix + XoY^ix + Ai)Pi ■ ■ ■ (x + A,r 

(where the Aj G F are mutually distinct and, for z = 0, G Z>o) we have 

the following decomposition of the unit in ¥[x]. We write Pi[x] := (x + Aj)^' and 
then P*[x] for the polynomial P[x]/ Pi, i = 0,1,- ■■ ,i. 

Lemma 3.1. There exist polynomials Qi[x], each of degree at most {pi — 1), so 
that 

1 = Qo[x]P'^[x] + Qi[x]P^[x] + ■■■ + Qe[x]P'[x]. 

Note that it is also easy to give a short inductive proof of this. The key specialisa- 
tion here is the bound on the degree of the QiS, otherwise the display is immediate 
from the polynomial variant of Proposition 12.71 

From this Lemma we immediately have specialisations of Theorem 12.21 and 
Corollary 12.31 However before we write these we would like explicitly to give 
formulae for the Qi[x] in the Lemma. We derive these in way which is rather 
suitable to our proposed applications. First observe that if B is an operator on V 
then, for p > 2, the solution space in V of B'^u = includes, for example, u such 
that Bu = 0. The solution space of B^u = is filtered. Given a solution u we may 
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obviously write m as a sum 

(27) u = + m(^) + ■ ■ ■ + 

where = 0, but such expansions are not unique. For example for any 

a G F we may take = (m — aBu) and u^^^ = aBu. 

Next observe that we may think of (V+\) as a nilpotent operator on the solution 
space Va of ([2]). Thus for A 7^ G F the operator (V + ii) is polynomially invertible 
on this space. With 

(28) A)-. 
we have 

(v + fir^\v + = tdv,. 

We shall write {V + /x)^^ to mean {V + fi)^^ o (V + fi)^^ and so forth. 

Now we construct the polynomial analogues of (l28l) . Write Q[x] to denote the 
polynomial (x + A)^. Suppose we consider F[x]/((5[a;]) meaning the algebra of 
polynomials modulo the ideal generated by Q[x]. As a multiplication operator on 
¥[x]/ {Q[x]), (x + A) is nilpotent and if, with fi as above, we write 

(29) (. + := - A)-(l + Ix^i + ■ ■ ■ + J^T^) 

then, 

(x + fi)x^{x + /i) = 1 mod {Q[x]). 
Similarly considering P[x] as in (l26ll note that in F[x]/(P[x]) the polynomial 
(x + Aj) is nilpotent as a multiplication operator on [ni^j=o(^ + -^i)^^]' since 
{x + Xif^ [ Yli^-=oi^ + XjY'] =P[x]. Consider now the vector in F[x]/ (-P[x]) given 
by, 

(30) Pr,[x] := [ n (a: + A.O^f ] [ ff (x + A,)^'^] . 

We will also view this as a multiplication operator on F[x]/(P[x]). Now from 
Lemma \3A\ we have 1 = ^i=o^«W-^*N ^[x]/ {P[x]) , where it should be noted 
we use the same notation for the polynomials 1 and Qi[x] and so forth as well as 
for their image in F[x]/(P[x]). Applying Prj[x] to both sides of this identity we 
have 

Pri[x]l = Pr4x]Qi[x]P*[x] mod (P[x]) 

since if /c G {0, ■■■,£} is distinct from i then Pr^xjP^'fx] vanishes modulo (P[x]). 
But now note that (x + Aj) is nilpotent on (5j[x]P*[x], in F[x]/(P[x]), as (x + 
Ai)P"P^[x] = P[x]. Thus if i ^ j G {0, ■ ■ ■ ,£} then, for example, (x + Aj)^/(x + A^) 
acts as the identity on Hence 

Pr4x]Q4x]P'[x] = Qi[x]P'[x] mod (P[x]), 

and so 

Fii [x] = Pr, [x] 1 = Qi [x] P' [x] mod (P [x] ) . 
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Thus from the Lemma and these observations we have 

(31) l = ^Pr4x] mod (P[x]). 

Finally we note that we may normalise the formula for Prjx]. Each term (x + 
in the product 

is a sum of powers of (x + Aj). In [x + Aj)^^^ it is only necessary to keep these 
powers up to (x + Ai)^'"^ as, recall, + Aj)^* [ni^j=o(^ + -^i)^^] ~ P\^\- Similarly, 
since we are applying the result to [ni^j^=o(^ + '^if'\ ^"^^ calculating modulo 
(P[x]), we may then expand the product ni^j=o(^ + writing the result as a 

linear combination of powers of [x + Aj) but always keeping only powers [x + \if 
for q G Z>o such that q < Pi — I. Let us write 

for this normalised formula for Qi [x] . Thus we have 

PrfW =iV([ n (x + A,),-f])[ n (x + A,r] 

for the corresponding normalised formula for Prj[x], i = 0, ■ ■ ■ ,£. So we have 

e j=i k=e 

i=E^([ n i^+^jfxn)[ n i^+^>^r] ^od{p[x]). 

But now observe that the normalised formula in the display has degree at most 
p — 1 in X, where p denotes the degree of P[x]. Thus we have the following result. 

Theorem 3.2. In W[x] we have the identity 

(32) 1 = X^Prf N = n + ^^■)aT1) [ n + • 
We associate this to the polynomial ^2E). 

(It seems likely that this identity is known from the theory of partial fractions.) 
It follows that for P of the form ([3]) we have the decomposition identity ([6j) with 

■■= N{ m;U(1^ + and := ng=o(^ + A.)'^ as follows. 

Corollary 3.3. Let V be a vector space over a field F. Suppose that V : V V is 
a linear endomorphism and that Aq, Ai, ■ ■ ■ , A^ G F are mutually distinct. We have 
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the identity in End(V); 

£ I j=£ k=£ 

= 5: prf = iv( [ n + a.oaTi ) [ n + ^^y-] 

where for i = 0, pj G Z>i. 

We obtain immediately the following specialisations of the results from Section 

EH 

Theorem 3.4. Let V be a vector space over a field F. Suppose that P : V — > V 
is a linear endomorphism and consider P : V — > V given by ^ with Xq, . . . , 
mutually distinct. Then there is a canonical and unique direct sum decomposition 
of the the null space for P, 

(33) Vp = ©LoVa, , 

where, for each i in the sum, Vx- is the solution space for (V + Aj)^\ This is 
executed by a canonical decomposition of the identity on Vp 

t 

idvp = ^Froii 

i=0 

where Projj : Vp V\^, z = 0, ■ ■ ■ , are projections given by the formula 

j=t k=£ 

(34) Proj, := [ n + ^^•)?] [ H + ^'^)''] ' 

Note in the theorem we have used the fact that we may omit the normalisation of 
(l30l) . since P[D] annihilates u. For the inhomogeneous problems Pu = f: 

Theorem 3.5. Let P be as above. Let us fix f E V. There is a 1-1 relationship 
between solutions u G V of Pu = f and solutions (mq, ■ ■ ■ ,ue) G ©^^^V of the 
problem 

(35) {V + XoY'uo = /,■■■, + X.pu, = f. 

Writing Vp for the solution space of Pu = f and (for i = 0, ■ ■ ■ , £) v£ for the 
solution space of {V + Ai)^'n = /. The map F : Vp x^^qV/. is given by 

j=£ j=e 

^^(11 (^^ + A,)^^■■■, n (^ + ^.)'0 , 

0^j=0 £=^j=0 

with inverse B : x^^qV£ — > Vp given by 

K- ■■,«.) ^ E ^( n + 

i=0 ij^j=0 
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Remark 3.6. Tuning our earlier discussion to the current setting we could opt to 
expand each Ui G Va^ with respect to the canonical filtration, say Ui = uf^ + ■ — h 

M-^' ^\ Although such expansions are not unique, we note here that the explicit 
form of the projection Proj^ given in ( l34l l gives such an expansion determined 

canonically by P. The point is this. Let us fix j 7^ i and write V^^ for the subspace 

of elements vectors h in Va^ satisfying (V + XiY^'^h = 0. First ni^j=o('^ + ^jY^"^ 

is in Va, = V|°\ Thus from ((28]) it follows that 



has the form 

/,(0) + . . . + hin-i) 

where h^''^ G V^f . Now {V + Ai) : v[f V^'^^^ where we view s e Zp.. Thus 
subsequent applications of {V + Xk)xl 7^ ^ ^ J) preserve this form and 
yield, in the end, an expression 



Ui = ai(l + ai,i{V + Xi) + ... + ai,p,_i(I? + A^)*''"^ 



n + >^^y' 



U 



where a, and ctij are determined explicitly by this process, and in fact it is easily 
seen that 

An important (generic) case of Theorems 13.41 and 13.51 is when we have 

P=(P + Ao)(P + Ai)---(P + A,) 
with the Aj mutually distinct. Then the situation simplifies as follows. 

Proposition 3.7. Let P be as in ^ with po = pi = ■ ■ ■ = pe = 1. Then for 
i = 0,--- ,i 

This follows immediately from the discussion in the remark above or is easily 
verified directly. 

For a linear operator P : V ^ V let us say that /i G F is in the spectrum of P 
(/i G Spec P) if (P — /i) : V — > V is not invertible (since we are not assuming that 
V is a Banach space). Suppose that P — fj, has the form ni=o(-^ + '^«) ('^ith the 
Aj G F not necessarily distinct) for some linear operator D : V — > V. Then, since all 
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factors commute, (P — /i) is injective (surjective) if and only if each of the factors 
{D + Aj) is injective (resp. surjective). Thus if F is an algebraically closed field 
and P = P[D] is polynomial in V then the spectrum of P is obviously generated 
by the spectrum oiV] jj, E Spec P if and only if /i = P[A] where A G Spec V. 

From the Theorem ll.il the eigenspaces are determined by the generalised eigen- 
vectors of V. We assume P to be a linear endomorphism operator on V, a vector 
space over an algebraically closed field F, in the summary here. 

Corollary 3.8. Let P = P[T>] be polynomial in V. Then {lJ.,u) is an eigenvalue, 
eigenvector pair for P if and only if for some A; G {1, ■ ■ ■ , deg(P)} 

U = Ui + ■ ■ ■ + Uk , 7^ Mj, z = 1, ■ ■ ■ , fc, 

where, for each i G {1, ■ ■ ■ , k}, (T> — X^Y'Ui = and Aj is a multiplicity Pi solution 
of of the polynomial equation (P — jj,)[x] = . 

Of course one could study generalised eigenspaces for P in the same way. 

3.2. The real case. If we work over a field that is not algebraically closed then 
the situation, in general, is different from Theorem ll.il since the polynomial P[x] 
may not factorise fully. However all is not lost. We illustrate the situation in 
the case that F is R, the field of real numbers. This case can be dealt with via 
complexification. By viewing P[x] as a polynomial in C[x] from the fundamental 
theorem of algebra we obtain a factorisation. 

(36) P[x] = (nz'o'ix + x^Y'^) (n:::='o'(a; + f^mY^-ix + f^mY^-) 

Here the — Aj G M, i = 0, ■ ■ ■ , £i, are the mutually distinct real roots and besides 
these there also the pairs of complex conjugate roots —Hm, —'^m G C \ R, with the 
Km mutually distinct for m = 0, ■ ■ ■ ,£2- So from Theorem 13.21 we have 

(37) 1 = E W + E W + W ) ' 

i=0 m=0 

where we have made an obvious adaption of the notation. By inspecting the 
formula there (i.e. (l32ll ) we see that Pr^[a;] is real, for 2 = 0, ■ ■ ■ , ii, and so is each 
sum (Prf \x] + Prf ), m = 0, - ■ ■£2- We note also that in each (Pi^ \x] + Prf ) 
there is a common factor P*"™ := P[x]/{x + kY'"^{x + kY""^ ■ Thus, in summary, 
by combining conjugate factors we obtain a real identity in R[x] of the form 

(38) 1 = ( E ^^^0 + ( E (Qm^p'^n) ' 

1=0 m=0 

where the Qi and the Qmm are obtained explicitly from (1371) . and where each term 
in the sum has polynomial degree less than the degree of P[x\. Thus from the 
general results of Section 12.11 we obtain the following. 

Corollary 3.9 (A real version of Theorem 13. 5p . Let V be a real vector space 
and P = PlV] : V — V an operator polynomial in V : V ^ V. Assume the 
complexification of P factors as in where Aj G M and Kj = jjLj + ii/j G C \ R 
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with /ij, Uj G M, Uj 7^ and the Aj 's and Kj 's are mutually distinct. Then the null 
space Vp, for P, admits a canonical and unique direct sum decomposition 

(39) Vp = 0Va.©0V,„., 

i=0 fc=0 

where jor each k in the sum, Vfi^^Uk is the solution space for (P^ + 2yUjD + /i| + z/|)'"-'. 

Fixing / G V, there is a 1-1 relationship between solutions m G V of Pu = f and 
solutions (mq, . . . , M^^, Mg, . . . , M^^) of the problem 

{V + Ao)^°«o = /,..., (I^ + \,T'u,, = f, 
{V' + 2^oV + + = f,...,{V' + 2fif,,V + fil + vlr^u',^ = f. 

The mappings relating these are given by F and B in Theorem\KM using ^EE). 

3.3. Operators polynomial in commuting endomorphisms. We now move 
to the general situation for this section. As above let us write V to denote a 
vector space over some field F. Suppose that : V — V, i = 1, ■ ■ ■ , fc, are 
non-trivial linear endomorphisms that are mutually commuting: T>iDj = VjVi for 
i,j G {!,■■■ ,k} . We obtain a commutative algebra W[D] of consisting of those 
endomorphisms V V which may be given by expressions polynomial (with 
coefficients in F) in the Vi. We write x = {xi, . . . ,Xk) for the multivariable 
indeterminate, and W[x] for the algebra of polynomials in the variables Xi, . . . ,Xk 
over the field F. Generalising the case of single variable polynomials, there is a 
unital algebra epimorphism from F[a:;] onto F[T>] given by formally replacing each 
variable Xj, in a polynomial, with Vi. 

Given polynomials Pq [^] > -Pi [^c] , ■ ■ ■ ,-P£[a;] G ¥[x] consider the product polyno- 
mial 

(40) P[x] = Po[x]P,[x]---P,[x]. 

With L = {0, 1, ■ ■ ■ ,L}, we carry over, in an obvious way, the labelling from 
Sections 12.21 and [2731 via elements of the power set 2^; products of the polynomial 
Pi[x] are labelled by the corresponding subset of L. For example for J C L, Pj[x] 
means rijeJ-^il*]' '^hile P"^[a;] mean Pl\j[x]. 

With a view to linking to the constructions above, we seek polynomials Qjfa;] G 
F[a3], J G a C 2^ satisfying the identity 

(41) i = J2QA^]p'[^] , 

or equivalently 

(42) 1 G {P-^[x] : J ea) 

where (..) denotes the ideal in ¥[x\ generated by the enclosed polynomials. Via 
the polynomial analogue of Proposition 12.71 we may equivalently study the "dual" 
problem of finding sets /5 C 2^ so that for each / G /3 we have 



(43) 



1 G {Pi[x\ : lel) 
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We may use algebraic geometry to shed light on this problem. Let us write 
Af{S[x]) for the algebraic variety determined by the polynomial S[x] G W[x] (i.e. 
{x G F^' I S[x] = 0}) and put A/} := Af{Pi) and Af^ := AA(P'^) for /, J G 2^. 
Clearly the condition (Il3l) requires Hi^jAfi = 0, because njg/A/"j C A/'(l) = 0. 
Comparing to (l43ll . the condition riig/A/i = is easier to verify, at least in simple 
cases, but it is generally weaker. However this depends on the field. In particular, 
it follows from the (weak form of) Hilbert's Nullstellensatz (see e.g. [5l Chapter 4, 
Theorem 1]) that if F is algebraically closed then 

IG {P^[x] : ^G/) ^ae,M = 0; 

i.e. for F algebraically closed the condition 1 G {Pi[x] : i G /) is equivalent to the 
polynomials Pi[x] : i G I having no common zero. (Note the previous display 
does not hold for F = R, e.g. take Pi[x,y] = + 1 and Pj[x,y] = y"^.) Thus with 
the notation introduced at the start of this section we have the following. 

Theorem 3.10. For P[X>] G W[D], (T> = (Pi,-- - ,Vk)) with F algebraically 
closed, 

(44) P[V] = Po[V]Pi[V]---Pe[D] 

is an algebraic dual P- decomposition o/ P[X'] if and only if, for the polynomials 
Pi[x] corresponding to the factors Pi['D], we have 

Hi^iAfi = 0, for all I e p . 

Generically for (3 G 2^, such that for all I & (3, \I\ > k, [44\ ) a dual (3- 
decomposition. 

Here, in an obvious way, we are using the term algebraic dual [3 -decomposition to 
mean a dual /^-decomposition that arises from the analogous polynomial identities 
as discussed. The last statement holds because generically Ht^jAfi has codimension 
|/|. 

The Theorem indicates immediately why one expects very strong results in the 
case of operators polynomial in a single operator. In one dimension algebraic 
varieties are generically disjoint. According to the theorem the situation is not 
much weaker for operators polynomial in several commuting operators. If we fix 
k then still we may say that generically operators polynomial in the operators 
Vq, ■ ■ - Vk admit (algebraic) a-decompositions. It is clear that Theorem 13.101 may 
be used to easily construct examples of all varieties of a-decompositions. 

We note that the results here are perhaps suggested by the general ideas of 
algebraic invertibility developed in |17| and references therein. However explicit 
links with the development in that source are currently far from clear. 

4. Symmetries 

Suppose that P is a linear endomorphism of a vector space V, over a field F. As 
above we write Vp for the kernel of P. Let us say that a linear map S" : V ^ V is a 
strong symmetry of P if S" preserves each of the eigenspaces of P. For example, if a 
V endomorphism S commutes with P, that is on V we have [S, P] := SP — PS = 0, 
then S* is a strong symmetry. On the other hand let us say that a linear operator 
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: Vp — i> V is a weak symmetry of P if S" has image in Vp C V. That is if S 
takes P-solutions to P-solutions. For example, ii S : Vp ^ V satisfies PS = S'P 
for some hnear operator S" : V ^ V, then S" is a weak symmetry. Evidently weak 
symmetries may be composed and via this operation yield an algebra. Similarly 
for strong symmetries. 

Given P as above, let us write Wp for the space of weak symmetries of P. In 
the case that P admits an algebraic decomposition P = PqPi ■ ■ ■ Pg (as in section 
12.11) then we obtain a corresponding decomposition of Wp, as a vector space. 
First one further item of notation. Let us write Wij for the vector space of linear 
homomorphisms H : Vj ^ Vj where, recall, Vj is the null space of Pj. Here we 
carry over notation from Section 12. 1[ 

Theorem 4.1. For P : V — > V, with ^ giving an algebraic decomposition, we 
have a canonical vector space decomposition, 



Proof. For H E Wij w obtain an element in Wp by forming H o Projj. This is 
inverted by the map taking arbitrary S G Wp to the composition 

Projj oS o Proj^- 



Note that for Hj^ G Wjk and Hij G Wij we have Hij o Hj^ G Wj^. Thus, 
identifying Wp with Q)Wij via the isomorphism in the Theorem, we see that for 
each i = 0,1, ■ ■ ■ ,£, Wu is a subalgebra of Wp. Evidently the algebra structure 
of Wp arises from that of these subalgebras plus the interlacing introduced by the 
spaces of homomorphisms W^, where i and j are distinct. Overall, understanding 
the algebraic structure of Wp is reduced to understanding the spaces Wij. 

Now suppose that F is an algebraically closed field and P is any polynomial in V. 
Recall from Corollary 13.81 that for a given G F the corresponding P-eigenspace 
(for simplicity of discussion we will allow this to be possibly trivial) V decomposes 

k 

into a direct sum = ©j^iVa, where the Va^ are generalised eigenspaces for V. 
Evidently we have the following observation. 

Proposition 4.2. If P : V V is a linear operator non-trivially polynomial in 
V and 5* : V — V preserves all generalised eigenspaces for V, then S is a strong 
symmetry for P. 

So for example any polynomial in V (viewed as a linear operator V ^ V) is a 
strong symmetry. 

The conditions in the Proposition are obviously too strict to generate all strong 
symmetries in general. It would be interesting to understand the precise relation- 
ship between strong symmetries for operators P, as in the Proposition, and the 
eigenspace information for T>. As a passing note we make a final observation in 
this direction. It is clear that if we fix /x in F then the restriction to of the 
linear maps S : V ^ V that preserve yields a space Wp_^ which is defined 
in the same way as the space of weak symmetries for the operator P — jj,, except 
that it consists of maps V ^ V (the domain is not taken to be the solution space). 



Wp = ©: 




in Wji. 



□ 
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Thus this may be analysed as was done for Wp above. The situation is rather sim- 
ple in lower degree cases. For example, the following proposition describes strong 
symmetries explicitly for P of degree 2. 

Proposition 4.3. Let P = {V + \i){T> + A2), with Ai, A2 G F (not necessarily 
distinct) and where T> -.V ^ V is a linear operator. For ^ G F, denote by V| the 
solution space of (T> + S,)^ for k E N. Then S : V V is a strong symmetry of P 
if and only if the following three conditions hold: 
(i) if —$,0 '■= — ^Ji^ G SpecD then S preserves V|j 

(a) if -(Ai + A2 - ^ Spec D, io then S preserves V| © '^x^+x^-^ 
(Hi) if G Spec D A — (Ai + A2 — Spec D then S preserves V|. 

Proof. Consider the decomposition of P — to irredicibles, i.e. 

P-/i= (P + Ai)(I) + A2)-/i= + + 
where /x G F and ^1, ^2 G F are not necessarily distinct. Then clearly ^1 + ^2 = 
A1-I-A2 and any pair ^1, ,^2 such that ^1+^2 = A1-I-A2 satisfies the previous display for 
some yU G F. Thus the strong symmetries are precisely linear mappings preserving 
the solution space of {V + ii){V + ^2) for every ,^1, ^2 such that ^1 + ^2 = Ai + A2. 
Using Theorem 13.71 the proposition follows. □ 

5. CONFORMAL LAPLACIAN OPERATORS AND EiNSTEIN MANIFOLDS 

On a smooth Riemannian or pseudo-Riemannian manifold {M"',g) let us write 
d for the exterior derivative and S for its formal adjoint (as in e.g. [I]). On the 
space of smooth /c-forms A'^ the form Laplacian is given by A = 6d + d6. Consider 
the operator 

Z = A - A^ 

where 7^ A G C. We may extend this to an operator on A*, the space of all 
smooth differential forms. Thus we have A = where D is the Dirac operator 
d + 6, hence 

Z= (D + \){D-X) 

on A*. Thus fixing / G A*, solutions m G A° of the problem Zu = f are in 1-1 
correspondence with solutions G (A*)^ of the problem 

{D + X)u+ = f {D-X)u^ = f, 

where we view A° C A*. The map from u, solving Zu = /, to a solution of the 
display is 

(45) {{D- A)m, {D + X)u) , 
while the inverse is 

(46) (m+,m_) ^ ^(m_ -M+). 

Using the grading of forms by degree, we may apply these tools to Z as an 
operator on functions A°. It is easily seen that (l46l l specialises to a map from 
A° (A° © A^)2, inverted by ([46]) as a map (A° © A^)^ A°. Fixing / G 
A°, this gives a 1-1 relationship between functions u solving Zu = f and pairs 
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e (A° © A^)^ solving {D + X)u+ = f and (D - A)m_ = /. In fact, once 
again using the grading of forms, one sees that the second order equation Zu = f 
is in fact equivalent to either one of the first order equations {D + \)u+ = f or 
(D - A)w_ = /. 

Operators of the form of Z arise naturally in Riemannian geometry. The con- 
formal Laplacian F : A° ^ A*^ is given by the formula 

An[n — 1) 

where Sc is the scalar curvature. Thus this is of the same form as Z on manifolds 
where Sc is constant and non-zero. 

A class of constant scalar curvature manifolds are the Einstein structures. A 
pseudo- Riemannian manifold (M, g) is said to be Einstein if its Ricci curvature 
is proportional to the metric (all structures will be taken to be smooth). We 
refer the reader to [I] for background on the meaning of these statements and the 
importance of Einstein structures. The conventions below follow [I2j except that 
we will use the "positive energy" Laplacian A as above (it may be also given as 
A = V*V, where V is the Levi-Civita connection and V* its formal adjoint. We 
assume the dimension of M to be at least 3. The GJMS conformal Laplacians 
of |15| are in general given by extremely complicated formulae, see [Hj. However 
on conformally Einstein manifolds we may choose an Einstein metric g. Then the 
order 2k GJMS operator may be viewed as an operator : —>■ A° and the 
formulae for these may be simplified dramatically. On Einstein n-manifolds the 
Pk is given by [l2] 



k 

(47) Pfc = J](A + qSc), 

i=l 

where Cj = (n + 2i — 2)(n — 22)/(4n(n — 1)) and Sc is the scalar curvature, that 
is the metric trace of the Ricci curvature. (For the standard sphere as a special 
case the formula (l47l ) was known to Branson [2].) On even manifolds the GJMS 
operators exist only up to order n. However for conformally Einstein structures 
it is shown in [l2] that, in a suitable sense, the family extends to all even orders. 
So for our current purposes for any k G Z>o we term the operator ( l47ll a GJMS 
operator. (We should also note that in line with our conventions for the sign of 
the Laplacian, the GJMS operator Pk as above is (— l)'^ times the corresponding 
operator in jlSj). 

Since the scalar curvature Sc is necessarily constant on Einstein manifolds it 
follows that Pk is polynomial in A and so we may immediately apply the results 
above to relate the null space of Pk with the generalised eigenvalues of the Lapla- 
cian. In the setting of compact manifolds of Riemannian signature it was noted in 
[T2] that we have such information via standard Hodge theory (or one could use 
functional calculus). The gain here is that we obtain related information in any 
signature and without any assumption of compactness. 
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The left (i.e. i = 1) factor in the expansion (l47j) is in fact the conformal Laplacian 
Y which plays a central role in spectral theory. So let us instead rephrase the 
Theorem 1.2 from [12] in terms of this. 



Theorem 5.1. On a pseudo-Riemannian n-manifold with Einstein metric, the 
order 2k GJMS operator is given by 

k 

(48) Pk = l[{Y + b,Sc), 

i=l 

where h = ^r—rr- 

' n(n—l) 

Note that when So 7^ the scalars 6jSc are mutually distinct. Thus from 
Theorem 13.41 and writing A/'(Pfc) for the null space of Pk as an operator on smooth 
real valued functions, we have the following. 

Theorem 5.2. On a pseudo-Riemannian Einstein n-manifold with So 7^ the 
null space of Pk has a direct sum decomposition 

M{Pu) = ©tiM(r) , 

where MiiY) is the eigenspace for Y with eigenvalue —6^80. 

Of course the machinery implies in the case of Sc = 0, but in this case the result is 
obvious: the null space is a generalised eigenspace for Y with generalised eigenvalue 
0, that is Af{Pk) = Af{Y^). In all cases the projection Af{Pk) Mi{Y) is given by 
(I3III . Similarly, the eigenspectrum of P^ is determined by Corollary 13.81 

Theorem 5.3. On a pseudo-Riemannian Einstein n-manifold, (/x, /) is an eigen- 
value, eigenfunction pair for the GJMS operator Pk if and only if for some m G 

/ = /i H h , 7^ /i, z = 1, ■ ■ ■ , m, 

where, for each i E {1, ■ ■ ■ ,ni}, (Y — XiY'fi = and Aj is a multiplicity Pi solu- 
tion of of the polynomial equation {Pk — = 0. (Here we consider Pk as the 
polynomial in Y , i.e. given by [4^)-) 

The inhomogeneous problems yield the obvious simplification to second order 
problems. 

Proposition 5.4. On a pseudo-Riemannian Einstein n-manifold, the inhomoge- 
neous problem PkU = f , for the GJMS operator Pk, is equivalent to the second 
order problem 

{Y + 6iSc)^ii = /,..., (r + bkSc)uk = f . 

From a solution {ui,--- ,Uk) of this problem we obtain, using bi = ^^j^^rj), the 
solution u of PkU = f as 

fn{n-l).k-is^rY^ 1 n 

""'^^^ gin 

In fact in odd dimensions and also in even dimensions n for the operators Pk<n/2 
we may further reduce to first order problems using the ideas at the start of this 
section. Via different Dirac operators there are variations on this outcome. 
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5.1. DifFerential Weak symmetries. It is clear that in any special setting the 
general idea of symmetries may be tuned somewhat. In particular, we shall do 
this for differential operators on pseudo-Riemannian manifolds. Suppose that now 
V is a space of smooth sections of some vector bundle over a pseudo-Riemannian 
manifold and P : V — > V is a differential operator. Then we shall say that a weak 
symmetry 5* of the differential operator P is differential if is given by a differential 
operator on V. That is 5* is differential weak symmetry of P means that it is 
a differential operator S : V ^ V such that it preserves the solution space of 
P. (This is slightly different from Section [4] where we defined weak symmetries 
only on the solution space of P.) Since the composition of differential operators 
yields a differential operator the differential weak symmetries form a subalgebra of 
the weak symmetries for P. Similar ideas apply to strong symmetries which may 
also be required to be differential. The key point is that provided the projection 
operators Projj (from Corollary 12. 3p are differential then the general results from 
section [4] carry over functorially to this category. 

In particular we illustrate this in the setting as above. Here we take V to be 
the space of smooth functions £ on an Einstein manifold M (of dimension at least 
3). Let us write W^-" for the space of linear differential operators S : £ ^ £ 
with the property that, upon restriction to Mj, S takes values in Afi{Y), that is 
S : Afj ^ J\fi{Y). The differential operators in S map between eigenspaces of 
the conformal Laplacian Y. From Theorem 13.41 and Theorem 14.11 we deduce the 
following. 

Theorem 5.5. On a pseudo-Riemannian Einstein n-manifold with Sc ^ 0, the 
space Wpj. of differential weak symmetries of the order 2k GJMS operator Pk has 
a canonical vector space decomposition, 

An obvious specialisation is to consider conformally flat spaces and locally (i.e. 
on a contractible manifold). Since the GJMS operators are conformally invariant, 
their solution spaces are conformally stable and one may study these by choosing 
a conformal scale that is congenial for the problem. For a current purposes a scale 
that achieves a constant non-zero curvature is ideal since then (on such Einstein 
structures) Theorem 15.21 applies. In particular we may apply Theorem 15.51 to 
study this conformal problem. In the setting of Euclidean space, Eastwood and 
Eastwood-Leistner P, [9] have studied the "higher symmetries" of the Laplacian 
and its square. These are differential weak symmetries S with the property that 
(in a choice of conformal scale) PS = S'P where S' : £ ^ £ is a differential 
operator. In this flat setting the Laplacian agrees with the Yamabe operator 
while the square of the Laplacian is the order 4 GJMS operator (which is usually 
termed the Paneitz operator). Since their theory is essentially conformal it should 
be an interesting direction to carry their results for the square of the Laplacian, 
in [9], onto a constant curvature conformally flat space and then relate these to 
our observations above. Our tools above provide an alternative approach to such 
higher order problems and also provide a route for studying the related questions 
on general conformally Einstein manifolds. 
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